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We investigate multi-wave mixing associated with the strongly pump depleted regime of induced 
modulation instability (MI) in optical fibers. For a complete transfer of pump power into the 
sideband modes, we theoretically and experimentally demonstrate that it is necessary to use a 
much lower seeding modulation frequency than the peak MI gain value. Our analysis shows that a 
record 95 % of the input pump power is frequency converted into the comb of sidebands, in good 
quantitative agreement with analytical predictions based on the simplest exact breather solution of 
the nonlinear Schrodinger equation. 


I. INTRODUCTION 

As well known, MI in a nonlinear dispersive medium 
described by the nonlinear Schrodinger equation (NLSE) 
leads to the exponential growth of perturbations, at the 
expense of an intense continuous wave (CW) pump back¬ 
ground [IH3- Since the very beginning of nonlinear wave 
propagation studies (see Ref. 0 for a review), it has 
been realized that MI (known also under different names, 
e.g., Benjamin-Feir instability for deep water waves i) 
is a universal pattern generation mechanism in different 
physical contexts such as plasma physics [l|, fluidody- 
namics 0 , and optics Despite fifty years of studies, 
the investigation of MI is still an extremely active field 
of interdisciplinary research. 

Initial MI studies essentially aimed at the observation 
of the early stage of exponential growth of the wave per¬ 
turbation spectrum. It is only recently that the main 
focus of the experiments has been moved to the fully 
nonlinear (or long-term) evolution of MI [341^ . In spite 
of the valuable pioneering approaches [3 t 9! such prob¬ 
lem remains open even theoretically, at least in its most 
general formulation, thus continuing to require consider¬ 
able efforts [isl - l^ . Indeed, knowing the long-term evo¬ 
lution of MI is a key issue for the understanding of many 
complex phenomena in p hysics, such as Fermi-Pasta- 
Ulam (FPU) recurrence [ol- flllfldl . higher order pulse- 
splitting [3, HH [13, homoclinic structures 0, Hj, [23| - 
rogue wave formation [25l - [27l| . the development of tur¬ 
bulence [13, and supercontinuum generation [13 . Since 
optical pulse propagation in fibers is well described by the 
NLSF, optical fibers provide a formidable test bed for the 
experimental investigation of all of these processes. 

Whenever MI is induced in the anomalous dispersion 
regime of a fiber by means of a relatively weak time- 
periodic perturbation with period 27r/D, the nonlinear 
evolution of MI leads, via a cascade four-wave mixing 
(FWM) process, to the appearance of a comb of harmon¬ 
ics at frequencies cuq ± uD, n = 1,2,3,... (cug indicates 
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pump frequency). Comb components exhibit a mono¬ 
tonic growth at the expense of the pump, but only up 
to a characteristic distance, say, Z^: at this point, the 
pump is maximally depleted. Further on, the power flow 
is reversed from the sidebands back into the pump, and 
subsequent cycles of periodic power exchange among the 
comb waves (or FPU recurrence) are typically observed. 

For a given fiber, the actual degree of maximum pump 
power depletion is strongly affected by the value of the 
initial modulation frequency D, as well as by the pump 
power P. The initial rate of frequency conversion (or 
MI gain) is well known to peak at a certain frequency, 
say ilpm- At this frequency, nonlinear phase matching 
occurs, i.e. the pump power induced phase shift cancels 
the linear dispersive mismatch: |/32|Dp^ = 2"fP (/?2 < 0 
and 7 are the fiber dispersion and its nonlinear coeffi¬ 
cient, respectively). 

However, setting the initial modulation frequency 
equal to flpm cannot guarantee that optimal frequency 
conversion from the pump also occurs in a regime where 
the pump becomes substantially depleted. Indeed, as the 
pump is depleted, the modulation frequency that satis¬ 
fies the nonlinear phase matching condition progressively 
shifts to lower values. The latter argument suggests that 
the maximum pump depletion may be increased, by set¬ 
ting the initial modulation frequency to a value which 
is lower than flpm- In such situation, instead of imme¬ 
diately mis-matching the mixing process, pump deple¬ 
tion may actually usefully drive the pump and sidebands 
towards phase-matching, thus effectively increasing the 
overall conversion efficiency [3l| . Despite the general fun¬ 
damental and applicative interest of wave propagation 
problems described by the NLSE, to our knowledge, the 
conditions that lead to the optimum transfer of energy 
from the pump to the sidebands have not been properly 
clarified, neither theoretically nor experimentally. 

In this paper we show that this question can find a sim¬ 
ple answer in terms of the so-called Akhmediev breather 
(AB) solutions of the NLSE 0 , i, [H [H, E, ■ More¬ 
over, we present a conclusive experimental evidence of 
this theoretical finding. As a matter of fact, ABs provide 
a full analytical description of the frequency comb gener¬ 
ation process in fibers. As we will show below, cascaded 
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FWM is responsible for a quantitatively important dis¬ 
crepancy with respect to a simple analytical prediction 
which may result from the truncated three-wave mix¬ 
ing (TWM) approach [23, HIl- Note that the 3-mode 
truncation was previously used for the modeling of fre¬ 
quency conversion experiments in fibers , carried out in 
the strongly pump depleted regime [H, • 

In order to validate the AB approach, in this work 
we report a careful characterization of the output FWM 
efhciency as the modulation frequency is varied. We 
also perform cut-back measurements, in order to portray 
the axial evolution of the pump and sidebands. This 
allows us to demonstrate that a record pump depletion 
of 95 % can be achieved, provided that MI is induced by 
a modulation frequency which is as much as 30% lower 
than the phase-matching prediction for an undepleted 
pump, in excellent quantitative agreement with the AB 
solution. 


II. THEORY 


Let us consider first the optimal frequency conversion 
problem from a theoretical point of view. Neglecting fiber 
loss, the dynamics of the depleted stage of MI in our 
experiments is well described by the NLSE 




( 1 ) 


where /32 is group velocity dispersion and 7 is the non¬ 
linear fiber coefficient. MI is induced by adding a weak 
in-phase modulation to the CW input pump in Eq. o 


zi(2 = 0,f) = yP -b V^exp(iwt)]. (2) 


Here P is the total power, t/q and rjg = I — r]Q are the 
pump and the signal input power fractions, and w is the 
pump-signal frequency detuning. The central question 
that we want to address is the following: what is the seed 
modulation frequency lu that guarantees the strongest 
(optimal) coupling of the input pump power into the 
signal-idler harmonic pair? In order to express this con¬ 
dition in universal form, it is convenient to employ di¬ 
mensionless units, i.e., we assume in Eq. O ^2 = -1, 
7 = I, P = I. This means that distance z, time t, and 
the field u are measured in units of the nonlinear length 
Zni = ( 7 ^*)”^, of the characteristic time Tq = \/|/32|^ni, 
and of -v/P, respectively. In this way, the solution of Eqs. 

depends on a single parameter only, namely, the 
dimensionless seed frequency w = flTo = Q^y\f32\/{"fP)- 
It is well known that, in these units, the linear stabil¬ 
ity analysis (LSA) of the CW solution of the NLSE ([1]) 
yields MI gain for 0 < to < 2, with a peak gain at 
uj = uipm, ojpm = flpmPo = being the normalized 
phase-matching frequency. 

A simple answer to the problem of the optimal con¬ 
dition for frequency conversion from the pump to the 
sidebands may be obtained by using the truncated (but 


fully nonlinear) TWM equations for the n = 0, ±1 
modes. By substituting the ansatz u{z,t) = Uq{z) + 
Mi(z) exp(— iwt) + u-i{z) exp{iujt) in Eq. ([T]), and ne¬ 
glecting higher-order sideband generation, leads to a self- 
contained set of ordinary differential equations (ODEs) 
for the pump and the first-order sideband amplitudes 
uo{z),u±i{z) [13, mi. As shown in Ref. [3l|, such equa¬ 
tions can be further reduced, in terms of the variables 
rj = I Mop and (j) = (j>i + (p-i — 2(j)o, to a one-dimensional 
integrable nonlinear oscillator. The equivalent particle 
motion is described in the phase plane by the ODEs 
drj/dz = dH/dcj), d(j)/dz = —dH/drj, with the Hamil¬ 
tonian 

H = 2r]\/ (I — r])^ — cos (j) + (w^ -b l)r] — 3?7^/2. (3) 

Here a = |m_i p— |mi p is a Manley-Rowe conserved quan¬ 
tity, which reduces to a = for the initial condition Eq. 
©• This model predicts that full conversion from the 
pump to the sidebands occurs whenever the phase plane 
points corresponding to the the initial condition and to a 
vanishing pump fraction rj = 0 belong to the same level 
curve of the Hamiltonian, or H{r]Q = 1 — rjs) = H{rj = 0) 
in Eq. This condition allows us to find a remarkably 
simple expression for the optimum input modulation fre¬ 
quency (for any given value of the initial power fraction 
of the signal rjs) 


( 4 ) 

Equation ([4]) predicts that, in the limit case of a small 
input signal (rjs "Cl), the modulation frequency for max¬ 
imum pump depletion is a factor two lower than the 
value obtained from the LSA of the undepleted pump, 
Wpm = \/2. Moreover, under the 3-mode approxima¬ 
tion, complete conversion is only asymptotically reached 
(that is, after an indefinitely long propagation distance), 
since t he p oint ry = 0 is a saddle point of the Hamil¬ 
tonian [3ll |. This situation is summarized in Fig. [ija): 
here we display the minimal pump fraction that is ob¬ 
tained at the normalized distance Zd = ZdjZni [i.e., 
r]min = vi^d)] when using the truncated TWM model. 
The input power fraction of the signal is equal to 3 % 
{r]s = 0 . 03,770 = 0.97). As it can be seen, full pump de¬ 
pletion occurs at the frequency ^twm predicted by Eq. 
(|4|): the corresponding maximum depletion distance Zd 
diverges to infinity as is approached. 

Equation ([4]) is qualitatively correct in predicting that, 
for obtaining maximum pump depletion, it is necessary 
to use a modulation uj which is substantially lower than 
the peak LSA gain value. However, this prediction is 
not quite in quantitative agreement with the NLSE so¬ 
lutions, owing to the presence of higher-order sideband 
pairs (n > 2). The coupling of power to the higher har¬ 
monics of the initial modulation becomes indeed stronger 
at modulation frequencies lower than phase-matching fre¬ 
quency, thus affecting also the flow of power towards the 
primary modes (n = ±1). 


Opt _ 

^TWM — 
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FIG. 1. (a) Predictions of the TWM truncated model: Dis¬ 
tance Zd and residual pump power fraction rj{zd) correspond¬ 
ing to maximally depleted pump versus modulation frequency 
uj, for a fixed input signal fraction of 3 % {rjs = 0.03). The 
dashed vertical line stands for from Eq. Q. (b,c,d) 

Full NLSE computation: (b) Evolution of Fourier modal frac¬ 
tion of power at frequency u) = = 1; (c) Residual pump 

power fraction vs. cu: AB (Eq. ([6]), solid line) compared with 
NLSE simulations (circles and crosses); (d) Corresponding 
distance Zd from NLSE simulations compared with approxi¬ 
mation (dashed lines) from Eq. ([5|). 


The exact frequency comb dynamics is generally de¬ 
scribed by doubly periodic (in t and z) solutions of the 
NLSE. However it has been shown that, for 1 < a; < 2, 
these general solutions remain sufficiently close to the 
AB solution d, d| (see also Refs. [H, [23 for a deriva¬ 
tion following a different method). The latter is the sim¬ 
plest full solution of the NLSE which is homoclinic to 
the background, i.e., its associated phase plane trajec¬ 
tory connects the background to itself after a full cycle 
of evolution (strictly speaking, the AB is a heteroclinic 
solution, since the trajectory connects two phase shifted 
background solutions). As it occurs in the TWM model, 
this cycle of evolution has an infinite period in the lon¬ 
gitudinal coordinate z. At the apex of the cycle, the AB 
describes a fully developed train of pulses, which corre¬ 
sponds to a maximally depleted background. 

Although the spectrum of the AB solution is symmet¬ 
ric around the pump, recent studies [H, [H, [l^ [s^ have 
shown that ABs may be used to approximate with rea¬ 
sonably good accuracy the nonlinear stage of induced MI 
with an asymmetric modulation seed, i.e., when a single 
sideband input condition replaces Eq. ©. This typically 
requires a small input signal fraction ? 7 s <C 1, and rela¬ 
tively large modulation frequencies (i.e., 1 < a; < 2 , so 
that no harmonics fall under the MI gain bandwidth). 

Quite interestingly, we could derive from the AB solu¬ 
tions a simple analytical condition for the optimum input 
modulation frequency, that leads to maximum pump de¬ 


pletion. Let us expand the AB at its apex (corresponding 
to the distance Zd) in a Fourier series 


u 


peak 

AB 


(t) 


{ui‘^/2 — 1 ) + y^l — a;2/4cos(a;t) 
a/ 1 — a;2/4cos(a;t) — 1 

n 


( 5 ) 


Here are the Fourier coefficients, which can be explic¬ 
itly calculated. We obtain 


|-peafe|2 




( 6 ) 


for the pump (n = 0 ), whereas for sideband modes ±n. 


\~peak\2^^2 


2 — oj 
2 oj 
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Equation ([^ implies that the pump is totally depleted 
at = 1. This is confirmed by the numerical sim¬ 

ulation of the NLSE o with initial condition © : see 
Fig. mb), where we report the evolution of the power 
fraction of the pump and the first four sideband pairs. 
Solving the NLSE at different modulation frequencies w 
confirms that the parabolic law of Eq. ([^ indeed pro¬ 
vides a quantitatively accurate description of the max¬ 
imally depleted pump in the whole range 1 < w < 2 , 
regardless of the initial power fraction of the signal. This 
agreement is displayed in Fig. [T](c)i where we compare 
the results of NLSE simulations (with two different in¬ 
put signal fractions ijs = 2 x 10“^, 300 x 10“^, or pump 
fractions rjo = 0.9998 and 770 = 0.97), to the analytical ex¬ 
pression [Eq. ([ 6 ])] . As it can be seen, slight discrepancies 
only appear for significantly high input signal fractions 
(see crosses for r]s = 0.03), and in the range of modula¬ 
tion frequencies well below the optimum value = 1 . 

Note also that, although the input signal fraction does 
not significantly affect the amount of maximum pump 
depletion (at least for 1 < a; < 2 ), 77 ^ strongly affects the 
maximum pump depletion distance Zd- The AB solutions 
also provide a reasonably raod estimate for Zd in terms 
of the following formula (30| (an alternative formula is 
reported in Ref. M) 


Zd = 




In 


/a;(l — a;2/4) 


( 8 ) 


Fig. [T](d) shows that Eq. ([5]) is indeed in good agreement 
with the numerical solution of Eqs. (fTHa . 

In summary, based on the AB solutions, one predicts 
that total pump depletion occurs precisely at = 1 . 
Such modulation frequency is higher than the prediction 
of the 3-mode truncation (i.e., (jJtwm — l/\/2), but it is 
still substantially lower (by a factor v^) than the phase¬ 
matching frequency Wpm = \/2. 

There is another difference between the results of the 
3-mode truncation and the AB theory which is worth to 
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FIG. 2. Experimental set-up. Cl-PM, polarization main¬ 
taining coupler; EDFA, erbium-doped fiber amplifier; EOM, 
electro-optic modulator; C2, coupler; ISO, isolator; SMF, 
Corning 1.1 km single mode fiber; OSA, optical spectrum 
analyzer. 

emphasize. In the former case, the condition for optimal 
pump depletion necessarily coincides with the condition 
for maximum power in the sideband modes. Conversely, 
the AB solution implies that the optimum power con¬ 
version from the pump into, e.g., the primary sidebands 
occurs at the modulation frequency 10^=1 = — 1 ~ 

1.24 7 ^ ■ This frequency is obtained by maximiz¬ 

ing the power fraction \ui\'^ in Eq. © [ in general each 
sideband order n peaks at a different frequency, which 
can be easily calculated from Eq. ©]. This is because 
at = 1 higher-order sidebands concur to the full 

depletion of the pump [even if ^ 96% of the power is 
contained in the modes with n < 3, the fractional con¬ 
tent of n = 2 ,3 modes is far from being negligible, see 
also Eig. ©b)]. Viceversa, by slightly increasing the fre¬ 
quency above = 1, the AB spectrum narrows down 
(i.e. the sideband modes decay faster for increasing n). 
In this case, even if the pump is not depleted completely, 
the fraction of power in the n = 1 sideband pair may still 
reach its largest value. 


III. EXPERIMENTAL RESULTS 

In order to experimentally investigate the optimal con¬ 
ditions for pump depletion in the MI process, we em¬ 
ployed the set-up reported in Eig. [2l To allow for rel¬ 
atively large modulation frequencies, we induced MI by 
an asymmetric seed. A strong CW-pump (A = 1560 nm) 
and a weak tunable CW-probe were combined and in¬ 
tensity modulated to generate 12 ns square pulses at 5 
MHz repetition rate. For increasing the pump pulse peak 
power, we used an erbium-doped fiber amplifier (EDFA2) 
before launching the modulated pump into a 1.1 km long 
Corning SMF28 fiber. For such a relatively short fiber 
length, the presence of fiber loss can be safely neglected. 

In order to fully characterize the MI dynamics, we per¬ 
formed two different sets of experiments. In the first set 
we measured at the fiber output the residual pump and 
signal power fractions, as the seed frequency detuning 
(or modulation frequency) A/ = r2/27r was varied. We 
compared such results with the measured MI gain curve. 
In the second set of experiments, we performed cut-back 



FIG. 3. Symbols: experiments. MI gain (red triangles), resid¬ 
ual pump (blue squares) and signal (magenta circles) fractions 
versus frequency detuning A/ = Q./{2ti). Dashed curves: 
simulations. The same colour code has been used. 


measurements and recorded the full output spectra every 
20 m. 

Fig. [3] illustrates the dependence of the fraction of 
total output power in the pump and in the signal, re¬ 
spectively, vs. their frequency separation. By using an 
extremely weak signal, the (undepleted pump) MI gain 
curve (red triangles in Fig. [31) was also measured. The 
MI gain profile shows that the phase matching frequency 
is ^fpm — 80 GHz, in good agreement with the theo¬ 
retical prediction (~ 84 GHz). Next we kept the pump 
power [P = 2.83 W) unchanged, and we increased the 
signal power in order to investigate the MI process in 
the strongly depleted pump regime. In this case, the 
signal power was accurately adjusted, so that the the 
fiber length was kept nearly equal to the maximum de¬ 
pletion length Zd- As shown in Fig. [31 the data (blue 
squares) clearly indicate that a minimum residual pump 
power of 5 % (95 % depletion) was achieved at A/ = 55.7 
GHz. This frequency is in good quantitative agreement 
with the estimate = 1 , which corresponds, in real 
world units, to A/ = Afpm/V^ = 59.4 GHz. The 
peak of the output signal fraction (magenta circles) was 
observed at A/ = 67.8 GHz, a value which agrees re¬ 
markably well with the predicted value — 1-24, or 
A/ = A/pm/1-24 = 67.7 GHz. For the sake of clarity, all 
of our results are summarized in Table 1. 


TABLE I. Results summary 


r 

Perfect 

phase 

matching 

Optimal pump 
conversion frequency 

Optimal signal 
conversion 
frequency 

TWM 

AB 

Normalized units 

LJPM = \/2 

Opt ^ 

^TWM - -J= 

Opt -1 

“ab = 1 

uj°Z{ -1.24 

Physical 

Units 

Measured 

(GHz) 

80 

55.7 

67.8 

Calculated 

(GHz) 

84 

= 42 

A/pm/V^ 
= 59.4 

67.7 


Gain (dB) 
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FIG. 4. Observed spectra versus wavelength and distance, 
as reconstructed from cut-back measurements: (a) / = 80.3 
GHz; (b) / = 55.7 GHz. 


In order to fully validate our measurement results, 
we compared them with numerical solutions of the 
generalized NLSE (GNLSE, which is typically used 
to simulate broaband pulse propagation and super¬ 
continuum generation [29j), using the available fiber 
data {^2 = 24.3 ps^/km, third-order dispersion /Ja = 
0.14 ps^/km,'j = 1.2W~^km~^, and the loss coefficient 
a = 0.2 dB/km) plus a CW pump and seed pair. The 
outcome of our simulations is superimposed (blue and 
magenta dashed lines) to the data in Eig. [3l As can be 
seen, an excellent agreement is obtained for the value of 
the optimal maximum depletion frequency. 

However, at variance with the experimentally observed 
minimal residual pump of 5%, simulations predict that 
full depletion should be observed (see also Fig. [J). Ad¬ 
ditional numerical simulations matching the pulsed (as 
opposed to CW) input pump show that the small resid¬ 
ual pump fraction is due to the finite roll-off factor of the 
non ideal square pump pulses. Indeed, averaging over the 
pump power profile is well known to lead to incomplete 
frequency conversion in the MI process, as previously re¬ 
ported in, e.g.. Ref. (l^ . 

In our second set of experiments, we recorded the evo¬ 
lution along the fiber length of the pump and signal power 
by means of a cut back experiment in steps of 20 m. 
These results are displayed in Fig. [4] (log vertical scale) 
for two specific values of the pump-signal frequency shift. 
The first value is close to the perfect phase matching 
frequency (A/ = 80.3 GHz, Fig. IHa)), while the sec¬ 
ond value is close to the optimal conversion frequency 
(A/ = 55.7 GHz, Fig. HKb)). As we can see, higher- 



Fiber length (m) 


FIG. 5. Measured pump and signal fractions against dis¬ 
tance obtained from cut-back measurements, for frequencies 
/ = 80.3 GHz (close to peak linear gain) and / = 55.7 GHz 
(optimum frequency for depletion), respectively. 


order sidebands, that start to grow at longer distances 
when compared with the primary sidebands (n = I), play 
a non-negligible role around the point of maximal deple¬ 
tion. In addition. Fig. [4Kb) shows that pump depletion 
is more pronounced around 1000 m. 

To gain a clearer insight into the longitudinal variation 
of the pump and signal power fractions. Fig. |5| compares 
their variation for the two different modulation frequen¬ 
cies of Fig. m Fig. |5| shows that, for A/ = 80.3 GHz, a 
relatively large (i.e., 27%) residual pump power fraction 
is observed at Zd — 1000 m. Gonversely, data obtained 
at A/ = 55.7 GHz show that a marked enhancement of 
the maximal pump power depletion occurs at Zd — 980 
m, in good agreement with the theoretical prediction. 


IV. SUMMARY 

In conclusion, we reported a clear experimental evi¬ 
dence that optimal power transfer from a pump wave 
to its sideband modes, as described by the fully nonlin¬ 
ear stage of the MI process in optical fibers, occurs for 
an initial modulation frequency A/ = A/p^/v^ that is 
well below the nonlinear phase-matching frequency A/p^ 
where the MI gain peaks. This result is important as it 
establishes a fundamental property of the nonlinear evo¬ 
lution of modulation instabilities in nonlinear dispersive 
media. Our observations have also widespread applica¬ 
tions to optical frequency conversion devices. Indeed, 
as it was earlier pointed out for vector Mis in birefrin- 
gent fibers [s^, and confirmed experimentally 
the concept of the phase-matching of parametric mixing 
processes is of limited use in the strongly depleted pump 
regime, unless it is suitably extended by means of a fully 
nonlinear large-signal theory. 










6 


FUNDING INFORMATION 

This work was partly supported by the ANR FOPAFE 
(ANR-12-JS09-0005), TOPWAVE (ANR-13-JS04-0004) 
and NoAWE (ANR-14-ACHN-0014) projects, by the 
’’Fonds Europeen de Developpement Economique 


Regional”, the Labex CEMPI (ANR-ll-LABX-0007) 
and Equipex FLUX (ANR-ll-EQPX-0017) through the 
’’Programme Investissements d’Avenir”. Funding from 
Italian Ministry of Research (grant PRIN 2012BFNWZ2) 
is also gratefully acknowledged. 


[1] A. Vedenov and L. I. Rudakov, Sov. Phys. Dokl. 9, 1073 
(1965) [russian version, Dokl. Akad. Nauk SSSR 159, 767 
(1964)]. 

[2] V. I. Bespalov, and V. I. Talanov, “Filamentary struc¬ 
tures of light beams in nonlinear liquids,” JETP Lett. 3, 
307 (1966). 

[3] T. B. Benjamin, and J. E. Feir, “The disintegration of 
wavetrains on deep water. Part 1: Theory,” J. Fluid 
Mech. 27, 417-430 (1967). 

[4] H. C. Yuen and W. E. Ferguson, “Relationship between 
Benjamin-Feir instability and recurrence in the nonlinear 
Schrodinger equation”, Phys. Fluids 21, 1275 (1978). 

[5] N. N. Akhmediev, V.M. Eleonoskii, and N.E. Kulagin, 
“Generation of periodic trains of picosecond pulses in an 
optical fiber: exact solutions,” Sov. Phys. JETP 62, 894 
(1985) [Zh. Eksp. Teor. Fiz. 89, 1542 (1985)]. 

[6] N. N. Akhmediev, and V. I. Korneev, “Modulation insta¬ 
bility and periodic solutions of the nonlinear Schrodinger 
equation,” Theor. Math. Phys. 69, 1089-1093 (1987) 
]Teor. Mat. Fiz. 69, 189-194 (1986)]. 

[7] K. Tai, A. Hasegawa, and A. Tomita, “Observation of 
modulation instability in optical fibers,” Phys. Rev. Lett. 
56, 135-138 (1986). 

[8] V. E. Zakharov and L. A. Ostrovsky, “Modulation insta¬ 
bility: The beginning”, Physica D 238 540-548 (2009). 

[9] G. Van Simaeys, Ph. Emplit, and M. Haelterman, “Ex¬ 
perimental demonstration of the Fermi-Pasta-Ulam re¬ 
currence in a modulationally unstable optical wave,” 
Phys. Rev. Lett. 87, 033902 (2001). 

[10] G. Van Simaeys, Ph. Emplit, and M. Haelterman, “Ex¬ 
perimental study of the reversible behavior of modula- 
tional instability in optical hbers,” J. Opt. Soc. Am. B 
3, 477-486 (2002). 

[11] J. Beeckman, X. Hutsebaut, M. Haelterman, and K. 
Neyts, “Induced modulation instability and recurrence in 
nematic liquid crystals,” Opt. Express 18, 11185 (2007). 

[12] J. M. Dudley, G. Genty, F. Dias, B. Kibler, and 
N. Akhmediev, “Modulation instability, Akhmediev 
Breathers and continuous wave supercontinuum gener¬ 
ation”, Opt. Exp. 17, 21497-21508 (2009). 

[13] K. Hammani, B. Kibler, C. Finot, P. Morin, J. Fatome, J. 
M. Dudley, and G. Mihot, “Peregrine soliton generation 
and breakup in standard telecommunications hber”, Opt. 
Lett. 36, 112-115 (2011). 

[14] K. Hammani, B. Wetzel, B. Kibler, J. Fatome, G. Finot, 
G. Millot, N. Akhmediev, and J. M. Dudley, “Spec¬ 
tral dynamics of modulation instability described using 
Akhmediev breather theory”, Opt. Lett. 36, 2140-2143 
( 2011 ). 

[15] A. Bendahmane, A. Mussot, P. Szriftgiser, O. Zerkak, G. 
Genty, J. M. Dudley, and A. Kudlinski, “Experimental 
dynamics of Akhmediev breathers in a dispersion varying 
optical fiber”. Opt. Lett. 39, 4490-4493 (2011). 


[16] M. Erkintalo, K. Hammani, B. Kibler, C. Finot, N. 
Akhmediev, J. M. Dudley, and G. Genty, “Higher-Order 
Modulation Instability in Nonlinear Fiber Optics”, Phys. 
Rev. Lett. 107, 253901 (2011). 

[17] A. Mussot, A. Kudlinski, M. Droques, P. Szriftgiser, and 
N. Akhmediev, “Fermi-Pasta-Ulam recurrence in nonlin¬ 
ear fiber optics: the role of reversible and irreversible 
losses,” Phys. Rev. X 4, 011054 (2014). 

[18] V. E. Zakharov and A. A. Gelash, “Nonlinear Stage of 
Modulation Instability”, Phys. Rev. Lett. Ill, 054101 
(2013). 

[19] V. E. Zakharov and A. A. Gelash, “Superregular solitonic 
solutions: a novel scenario for the nonlinear stage of mod¬ 
ulation instability”. Nonlinearity 27 R1-R39 (2014). 

[20] G. Biondini and E. Fagerstrom, “The integrable nature 
of modulational instability”, SIAM J. Appl. Math. 75, 
136-163 (2015). 

[21] A. Calini and G. M. Schober, “Homoclinic chaos increases 
the likelihood of rogue wave formation”, Phys. Lett. A 
298, 335-349 (2000). 

[22] S. Wabnitz, and N. Akhmediev, “Efficient modulation 
frequency doubling by induced modulation instability”. 
Opt. Gommun. 283, 1152-1154 (2010). 

[23] M. J. Ablowitz, and B. M. Herbst, “On homoclinic struc¬ 
ture and numerically induced chaos for the nonlinear 
Schrodinger equation”, SIAM J. Appl. Math. 50, 339- 
351 (1990). 

[24] S. Trillo and S. Wabnitz, “Dynamics of the nonlinear 
modulational instability in optical fibers,” Opt. Lett. 16, 
986-988 (1991). 

[25] A. Osborne, “The random and deterministic dynamics of 
rogue waves in unidirectional, deep-water wave trains”, 
Marine structures 14, 275-293 (2001). 

[26] M. Onorato, S. Residori, U. Bortolozzo, A. Montina, and 
F. T. Arecchi, “Rogue waves and their generating mech¬ 
anisms in different physical contexts”, Phys. Rep. 528, 
47-89 (2013). 

[27] J. M. Dudley, F. Dias, M. Erkintalo, and G. Genty, “In¬ 
stabilities, breathers and rogue waves in optics”. Nature 
Photonics 8, 755-764 (2014). 

[28] M. J. Ablowitz, J. Hammack, D. Henderson, and G.M. 
Schober, “Modulated periodic Stokes waves in deep wa¬ 
ter,” Phys. Rev. Lett. 84, 887-890 (2000). 

[29] J. M. Dudley, G. Genty, and S. Coen, “Supercontinuum 
generation in photonic crystal fiber”. Rev. Mod. Phys. 
78, 1135-1184 (2006). 

[30] M. Erkintalo, G. Genty, B. Wetzel, and J. M. Dudley, 
“Akhmediev breather evolution in optical fiber for real¬ 
istic initial conditions”, Phys. Lett. A 375, 2029 (2011). 

[31] G. Cappellini, and S. Trillo, “Third-order three-wave 
mixing in single-mode fibers: exact solutions and spa¬ 
tial instability effects,” J. Opt. Soc. Am. B 8, 824-840 
(1991). 



7 


[32] M. E. Marhic, K. K. Y. Wong, M. C. Ho, and L. G. 
Kazovsky, “92% pump depletion in a continuous-wave 
one-pump fiber optical parametric amplifier”. Opt. Lett. 
26, 620-622 (2001). 

[33] P. Kylemark, H. Sunnerud, M. Karlsson, and P. A. An- 
drekson, “Semi-analytic saturation theory of fiber optical 
parametric amplifiers”, J. Lightwave Technol. 24, 3471- 
3479 (2006). 

[34] G. Gappellini and S. Trillo, “Energy conversion in de¬ 
generate four-photon mixing in birefringent fibers”. Opt. 
Lett. 16, 895-897 (1991). 


[35] S. Trillo, G. Millot, E. Seve, and S. Wabnitz, “Fail¬ 
ure of phase matching concept in large-signal paramet¬ 
ric frequency conversion”, Appl. Phys. Lett. 72, 150-152 
(1997). 

[36] G. Millot, E. Seve, S. Wabnitz, and S. Trillo, “Observa¬ 
tion of a novel large-signal four-photon instability in op¬ 
tical wave mixing”, Phys. Rev. Lett. 80, 504-507 (1998). 

[37] E. Seve, G. Millot, and S. Trillo, “Strong four-photon 
conversion regime of cross-phase modulation induced 
modulational instability”, Phys. Rev. E 61, 3139-3150 
( 2000 ). 


